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The description of excited state dynamics in energy transfer sys- 
tems constitutes a theoretical and experimental challenge In mod- 
ern chemical physics. A spectroscopic protocol which systemat- 
ically characterizes both coherent and dissipative processes of 
the probed chromophores is desired. Here, we show that a set 
of two-color photon-echo experiments performs quantum state 
tomography (QST) of the one-exclton manifold of a dimer by re- 
constructing its density matrix In real time. This possibility in 
turn allows for a complete description of excited state dynamics 
via quantum process tomography (QPT). Simulations of a noisy 
QPT experiment for an inhomogeneously broadened ensemble 
of model excltonic dimers show that the protocol distills rich In- 
formation about dissipative excltonic dynamics, which appears 
nontrlvially hidden in the signal monitored in single realizations 
of four-wave mixing experiments. 

nonlinear spectroscopy | quantum information processing | excitation energy 
transfer 

Abbreviations: QPT, quantum process tomography; QST, quantum state tomog- 
raphy; PE, photon echo 

Excitonic systems and the processes triggered upon their interac- 
tion with electromagnetic radiation are of fundamental physical 
and chemical interest 1 1 , 2, 3 , 4, 5 , 6|. In nonlinear optical spec- 
troscopy (NLOS), a series of ultrafast femtosecond pulses induces 
coherent vibrational and electronic dynamics in a molecule or nano- 
material, and the nonlinear polarization of the excitonic system is 
monitored both in the time and frequency domains (7l[8l. To interpret 
these experiments, theoretical modeling has proven essential, framed 
within the Liouville space formalism popularized by Mukamel |7|. 
Implicit in these calculations is the evolution of the quantum state of 
the dissipative system in the form of a density matrix. The detected 
polarization contains information of the time dependent density ma- 
trix of the system, although not in the most transparent way. An 
important problem is whether these experiments allow quantum state 
tomography (QST), that is, the determination of the density matrix 
of the probed system at different instants of time l9l 1101 . A more 
ambitious question is if a complete characterization of the quantum 
dynamics of the system can be performed via quantum process to- 
mography (QPT) 1 1 1 , 12 1, a protocol that we define in the next sec- 
tion. In this article, we show that both QST and QPT are possible for 
the single-exciton manifold of a coupled dimer of chromophores with 
a series of two-color photon echo (PE) experiments. We also present 
numerical simulations on a model system and show that robust QST 
and QPT is achievable even in the presence of experimental noise as 
well as inhomogeneous broadening. 

This article provides a conceptual presentation, and interested 
readers may find derivations and technical details in the Supporting 
Information (SI). 



the system is p{T) = Tr_B ptotai(r), where the trace is over the 
degrees of freedom of the bath. If the initial state is a product, 
Ptotai(O) — p(0) ® ps(0) (always with the same initial bath state 
ps (0)), then the evolution of the system may be expressed as a linear 
transformation LT3J : 

Pa6(r) = ^Xai,cd(r)pcd(0). [1] 

cd 

The central object of this article is the process matrix x{T), which is 
independent of the initial state p(0). As opposed to master equations 
which are written in differential form, Eq.[T]can be regarded as an in- 
tegrated equation of motion for every T. It holds both for Markovian 
and non-Markovian dynamics of the bath, and it always leads to pos- 
itive density matrices. Note that x{T) completely characterizes the 
dynamics of the system. Preserving Hermiticity, trace, and positivity 
of p(T) imply, respectively, the relations (SI-I) 

Xabcd{T) = Xbadc{T), [2] 

'^Xaacd(T) = Scd, [3] 
a 

^ Z*cXabcdiT)zbd > 0, [4] 
abed 

where z is any complex valued vector. Using Eqs.|2]and[3l for a sys- 
tem in a d-dimensional Hilbert space, x{T) is determined by — 
real valued parameters 1 1 1 1. Operationally, QPT can be defined as an 
experimental protocol to obtain x(T)- Spectroscopically, xiT) may 
be reconstructed by measuring p{T) (i.e., performing QST) given 
some choice of initial state p(0), where the p(0) are chosen succes- 
sively from a complete set of initial states 1741 1151 [T6l II7I . Since 
we are interested in energy transfer dynamics, this procedure shall 
be performed at several values of T. In this article, we show how to 
perform QPT for a model coupled excitonic dimer using two-color 
heterodyne photon-echo experiments. 



Description of the system and its interaction with light 

Consider an excitonic dimer interacting with a bath of phonons. The 
excitonic part of the Hamiltonian, describing the system with the en- 
vironment frozen in place, is given by l4l[8l lI8l : 
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Basic concepts of QPT 

Consider a quantum system which interacts with a bath. We can 
describe the full state of the system and the environment at time 
T by the density matrix ptotai(r). The reduced density matrix of 
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= iOaCtCa + iiJl3C'^Ci3. [5] 

where af and c"'' (a^ and Cj) are creation (annihilation) operators 
for site i £ {A, B} and delocalized j G {a, /?} excitons, respec- 
tively. UJA are the first and second site energies, J 7^ is the 
Coulombic coupling between the chromophores. We define the aver- 
age of site energies lo — \ (cja + i^b ) , difference A — i(ajA — i^s), 
and mixing angle Q = \ arctan (■^) • Then = cos SaA+sin Sas, 
C/3 = —waQaA + cos 9aB, Ua — uo + A sec 26, and oj^ = 
(jj — A sec 26. For convenience, we define the single exciton states 
\^) = CqIp), 1/3) = cjflg), where \g) is the ground state, and the 
biexcitonic state |/) = a\a'^\g) = CaC^\g)- The model Hamilto- 
nian does not account for exciton-exciton binding or repulsion terms, 
so the energy level of the biexciton is oj/ — oja + i^p = ^a + l^b- 



Denoting Ui 



, we have LOag = iOfp and iopg 



We are interested in the perturbation of the excitonic system due 
to three laser pulses: 



V{t') = -A^A- ei£(t'-t»){e''=- " 



-iuj-lit —t-i) 



+ C.C.}, [6] 



yields a signal that is proportional to the components of Pp_E(t') • 64 
oscillating at the frequencies of the LO, which is centered about 0)4, 
Q. In this excitonic model, the only optically allowed transitions are 
between states differing by one excitation, so the only nonzero matrix 
elements of fi are /x^j = fi^^ for ij — ag, fig, fa, fj3 (SI-II). 

In the following section, we present the main results of our study. 
We show that a carefully chosen set of two-color rephasing PE ex- 
periments can be used to perform a QPT of the first exciton mani- 
fold (Fig. 1). The preparation of initial states is achieved using the 
first two pulses at ti and t2- Initial states spanning the single exci- 
ton manifold are produced by using the four possible combinations 
of two different carrier frequencies for the first two pulses. In the 
terminology of PE experiments, these pulses define the so-called co- 
herence time interval r = t2 — ti. The time interval between the 
second and third pulses, called the waiting time T = ts — t2, defines 
the quantum channel 111] which we want to characterize by QPT. Fi- 
nally, we carry out QST of the output density matrix at the instant 
t-^. This task is indirectly performed by using the third pulse to se- 
lectively generate new dipole-active coherences, which are detected 
upon heterodyning with the fourth pulse at t4, that is, after the echo 
time t — t4 — has elapsed. Varying the third and fourth pulse fre- 
quencies yields sufficient linear equations for QST. This procedure 
naturally concludes the protocol of the desired QPT. 



where A is the intensity of the electric field, assumed weak, p, is the 
dipole operator, and ei,ti, ki,uji denote the polarization^ time cen- 
ter, wavevector, and carrier frequency of the i-th pulse. E{t') is the 
slowly varying pulse envelope, which we choose to be Gaussian with 

fixed width a for all pulses, E{t') = e^* ''^'^ . The polarization 
induced by the pulses on the molecule located at position r is given 
by ^(r', t') — Tr(/ip(r, t')). This quantity can be Fourier decom- 
posed along different wavevectors as P{r,t') = Ps{t')e^'°""^ , 
where the ks are linear combinations of wavevectors of the incoming 
fields. Radiation is produced due to the polarization (proportional to 
iP{r, t')). We can choose to study a single component Ps by detect- 
ing only the radiation moving in the direction fcs. This is achieved 
by interfering the radiation with a fourth pulse moving in the direc- 
tion ks, called the local oscillator (LO) |7|. In particular, we shall 
be interested in the time-integrated signal in the photon-echo (PE) 
direction, ki = kpE ~ — fei + fc2 + fe3. This heterodyne-detected 
signal [SpE]el',e2.'e3,'e!t , whcrc the subscripts indicate the polariza- 
tions of the four light pulses and the superscripts indicate their carrier 
frequencies, is proportional to 



/oo 
-00 



44)64 • iPpE{t') 



[7] 



Spatial integration over the volume of probed molecules selects out 
the component PpE{t') from the P{r,t'). The time integration 




(c) 



Results 

For purposes of the QPT protocol, we assume that the structural pa- 



rameters UJag, UJfig, H^g, Mag, My, 



Figure 1. A set of photon echo experiments can be regarded as a QPT. Pulses 
are centered about ti, (2, ts, t4. Time flows upwards in the diagram. The four 
pulses define the coherence (r), waiting (T) times, and echo [t) times. This ex- 
periment, in the language of quantum information processing, can be regarded as 
consisting of three stages: initial state preparation, free evolution, and detection 
of the output state of the waiting time. 
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are all known. Infor- 



mation about the transition frequencies can be obtained from a linear 
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Figure 2. Possible state preparations and detections. We list all the possible 
preparations and detections of elements of the density matrix at the waiting time 
T via a rephasing PE experiment. The double sided Feynman diagrams above list 
all the possible processes detected in a rephasing PE experiment. Each diagram 
is related to an element Xijqp{T) due to the prepared state \ q){p\ — 5pq\g){g\ 
at the beginning of the waiting time and the detected state ij at the end of it. By 
combining the preparations in (a), (b), (c), and (d), with the detections in (e), (f), 
(g), and (h), sixteen different types of processes can be enumerated, which can 
be classified according to the pulse frequencies of the associated perturbations. 



We use the word polarizoiion in two different ways: To denote (a) the orientation of oscillations 
of the electric field and (b) the density of electric dipole moments in a material. The meaning 
should be clear by the context. 

-More precisely, the monitored signal is proportional to /'^^ dfe*['^'i(* -^4 ) + ¥'] — 

44)64 ■ iP pE{t') + C.C., where two experiments are conducted by varying the phase ^ of 
the LO with respect to the emitted polarization to extract the real and imaginary terms of Eq. 
[71 For purposes of our discussion, it is enough to consider the complex valued signal. 
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absorption spectram, whereas the dipoles can be extracted from x- 
ray crystallography f 19 |. As shown in recent work of our group, with 
enough data from the PE experiments, it is also possible to extract 
these parameters self-consistently 1201 1211 . We proceed to describe 
the steps of the PE experiment on a coupled dimer that yield a QPT. 

Initial state preparation Before any electromagnetic perturbation, 
the excitonic system is in the ground state p(— oo) — \g){g\- Af- 
ter the first two pulses in the fei, fe2 directions act on the system, the 
effective density matrix p%l]t? (0) (at T = 0) is created. This density 
is second order in A and, combined with the third and fourth pulses, 
directly determines the signal. By applying second order perturbation 
theory and the rotating-wave approximation (RWA), we can define an 
effective initial state (Fig. 2a-d, SI-III): 

>^Qav{^){\q){p\-5p,\g){g\). [8] 
This state evolves during the waiting time T to give 

-pilTAT) = x(r)p^j:e"^(o), [9] 

which holds for T > 3a, that is, after the action of the first two pulses 
has effectively ended. Eq.|9]is of the form of Eq.[T] and therefore ap- 
pealing for our QPT purposes. The purely imaginary coefficients C£ . 
for p £ {a, /?} are proportional to the frequency components at cupg 
of the pulse which is centered at cji : 

Cl. = _AV2^e-"'("''«-"')'/', [10] 
and the propagator of the optical coherence is 

e.,(r) = e(r)e(-^"--^-)^ [11] 

which, for simplicity, has been taken to be the product of a coherent 
oscillatory term beating at a frequency lOij and an exponential decay 
with dephasing rate Fij, assumed to be known. 0(r) is the Heavi- 
side function, so the propagator is finite only for times r > 0. We 
have kept only the — fei and +fe2 components because those are the 
only contributions to the signal at fcp_B . 

Eq. ([8} has a simple interpretation, and can be easily read off 
from the Feynman diagrams depicted in Fig. 2a-d. Since we will be 
selecting only light in the kpE direction, we keep only the portion of 
the first pulse proportional to e*"i''^'i' and only the portion of the 
second pulse proportional to e^"^^''^'^' in Eq. |6l Then, since the 
state before any perturbation is \g){g\, the first pulse can only reso- 
nantly excite the bra in the RWA |7|, generating an optical coherence 
\g)kp\ (where p G {a, /?}) with amplitude CZ^ - This coherence un- 
dergoes free evolution for time r under Qgp(T) before the second 
pulse perturbs the system. In the RWA, this second pulse can act on 
the ket of \g) {p\ to yield \q) {p\ with amplitude and on the bra to 
create a hole ~\g){g\ with amplitude C£, , producing Eq.[8] The am- 
plitude of this prepared initial state is proportional to the alignment 
of the corresponding transition dipole moments with the polarization 
of the incoming fields. Once the initial state is prepared, it evolves 
via x{T), which is the object we want to characterize. A possible 
problem is the contamination of the initial states by terms propor- 
tional to a hole— every time there is a single-exciton manifold 
population \p){p\. This is not a difficulty if we assume: 

XabggiT) = SagStg, [12] 

that is, if the ground state population does not transform into any 
other state via free evolution. This is reasonable since we may ignore 
processes where phonons can induce upward optical transitions. We 
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similarly neglect spontaneous excitation from the single to double 
exciton states. 

One can see from Eq. [8] that a set of four linearly independent 
initial states can be generated by manipulating the frequency compo- 
nents of the pulses through C£j and C^., . It is sufficient to consider 
a pulse toolbox of two waveforms which create \a} and j/?) with dif- 
ferent amplitudes. For instance, by centering one waveform at uu+ in 
the vicinity of tj^g and the other at tj_ close to uj/jg , we can simul- 
taneously have: 

= C'f_ = c', 
C£+ = C"_ = C", [13] 

for purely imaginary numbers C' 7^ C" . The conceptually simplest 
choice, which we shall denote the maximum discrimination choice 
(MDC), and which is best for QPT purposes, is two waveforms each 
resonant with only one transition, so C' 3> C", and we can ne- 
glect C". Four linearly independent initial states can be prepared by 
choosing the waveform of each of the first two pulses from this pulse 
toolbox. 

Evolution The system evolves during the waiting time T after the ini- 
tial state is prepared. Transfers between coherences and populations 
are systematically described by x[T). The components of p%l]t2 
evolve in time, described by x{T), without assuming any particu- 
lar model for the bath or system-bath interaction, other than Eq. [12] 
By definition, the amplitude of the ij component of p^l'X^^iT) is 
proportional to Xiiqp{T), with the exception of the gg component, 
which is proportional to Xggpq{T) + Spq due to the contamination 
of the hole in the initial state. We note that de-excitation transfers 
from the single-exciton manifold to elements involving the ground 
state (gg, ga, gP) are expected to be small, since such processes are 
also unlikely to occur on femtosecond timescales due to either the 
phonon or the photon bath. Detailed analysis shows that our protocol 
can detect decay into gg but not into ga or gfi. We shall keep these 
Xaaqp{T) terms in our equations in order to monitor amplitude leak- 
age errors from the single-exciton manifold, providing a consistency 
check for treating the single-exciton manifold as an effective TLS in 
the timescale of interest. 

Detection The last two pulses provide an indirect QST of the state af- 
ter the waiting time. The third perturbation along -ffcs and centered 
at time ts will selectively probe certain coherences and populations 
of pelle2 (T)- As an illustration (see Fig. 2e,f), the component of the 
pulse matching the transition energy uj^g = tjJ/s and proportional 
to C^gWill, in the RWA, promote the resonant transitions \g) — >■ |a) 
and 1/3) — >■ I/) on the ket side, and the conjugate resonant transitions 
(q| — >■ {g\ and (/| — > on the bra side, the latter of which are 
irrelevant as we are ignoring transfers to the biexciton state during 
the waiting time. These transitions will generate two sets of opti- 
cally active coherences in the echo interval: ag, fj3 which oscillate 
at frequency ujag, and j3g, fa which oscillate at ujpg. These sets 
generate a polarization which interferes with the LO yielding sig- 
nals proportional to C"^ and C£^, respectively. The propagator for 
the echo time is taken to be as in Eq. [TT] A similar analysis can 
be repeated for the Cf^ term. Analogously to the preparation stage, 
the same toolbox of two different waveforms for the third and the 
fourth pulses allows discrimination of all final states of Pe^'eg^l^). 
Fig. 2 depicts double-sided Feynman diagrams for all possible com- 
binations of preparations and detections with four pulses each chosen 
from two waveforms, yielding sixteen experiments. By keeping track 
of these processes, the signal [SpE\'il]e2'ei^e!i rnay be compactly 
written as (SI-IV): 

[SpE]'il',e2'X3,'Bi (t, T, t) 

^ ^ Cf^-^C^^C^.^C^^P^'^^'g^^^^^g^{T,T,t), [14] 

p,q,r,s 
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where the proportionahty constant is purely real, and the expression 
holds for T,t> 3a. The terms Pii'e'atea.ei {r, T, t) are loosely po- 
larizations (in fact, they are proportional to i times polarizations) Q 
and are given by 



pp.q,a.a , j, 

ei ,62 ,e3 ,64 V ' 7 5 



^{Xsgqp{T) — Spq — XaaqpiT)) 

+ (A'//3 ■ e3)(Pfp ■ e4)gffi{t)xi3i3qp{T)]}, [15] 



and 



pP,q,a,^ ( y ^\ 
^ 61,62,63,64 V ' ' ' 

= ~(Pp3 ■ e^)itJ-qg ■ e2)Qqp{T) 

x{[((M//3 • e3)(My„ ■ e4)gfc{t) 

-iPcg ■ e3){fJ-i3g ■ ei)g,3g{t))xpc,qp(T)]}. [16] 

The remaining terms Pl[':£,^.„^,{T,T,t), Pl\''. t% {T,T,t) 
follow upon the interchange a -H- p. Eqs. 1 141161 are the 
main result of this article. Each Pe'^Je2%3,e4, represents the ob- 
served signal if the first (second,third,fourth) laser pulse is reso- 
nant only with the p (q,r,s) transition. The total measured signal 
[gpg]6^'6 f6"{.6"/ (r,r,t) is a weighted sum of these P^.^.-^iU,.,,. 
Equations Il5|16l show that each Pe'i.e2%3-e4 ^ linear combi- 
nation of elements of xiT), corresponding to the prepared initial 
states and detected final states. After collecting the sixteen signals 
[S'pEjej'jej, 63^,64'' with sach pulse carrier frequency a;,; chosen from 
{ijj+,uu-} as in Eq. J13b . with fixed polarizations e^, Eq. J 14b can 
be inverted to yield the elements of xiT) associated with the single- 
exciton manifold of the dimer, hence accomplishing the desired QSTs 
and QPT at once (see Table 1). Notice that in principle, for a given 
value of waiting time T, the one-dimensional (ID) measurements as- 
sociated with a single set of r, t values is enough for purposes of QPT 
of the single exciton manifold. In the most typical measurements, the 
sample has isotropically distributed chromophores, so Eq. 1 141 must 
be modified to include isotropic averaging (■)iso (SI- VI). Since the 
present QPT protocol does not rely on different polarization settings, 
we will assume for simplicity that each of the sixteen experiments 
is carried out with (ei, 62, 63, 64) — {z, z, z, z). Further technical 
details of the QPT protocol can be found in the next section as well 
as in SI-V-XII. 

Important observations (a) Difference between a standard PE ex- 
periment and QPT. In the current practice of NLOS, model Hamil- 
tonians with free parameters for the excitonic system, the bath, and 
the interaction between them are postulated, and the experimental 
spectra are fit to the model via calculation of response functions, 
from which structural and dynamical information is extracted |22 1. In 
our language, such experiments involve fitting models to complicated 
combinations of quantum processes associated with x{T)- QPT re- 
quires only a model for the excitonic system, but not for the bath 
or the system-bath coupling, making it suitable for probing systems 
where the bath dynamics are unknown. By definition, QPT extracts 
the elements of x{T), allowing a straightforward analysis of pro- 
cesses directly associated with the density matrix, such as dephasing 
and relaxation. 

(b) QPT can also be performed with control of time delays 
T, t instead of frequency control. Although ID measurements 
suffice for QPT, suppose the signal is collected for many values of 
r and t. Upon appropriately defined Fourier transformations of the 
signal along these variables, a two-dimensional electronic spectrum 
(2D-ES) can be constructed where the coherence propagators of Eq. 
dl lb manifest as four resonances about iVag and ujpg along both axes 



An important observation follows: the frequencies of the 
coherent evolutions in the coherence and echo times are the same as 
the frequencies of the first transition and the LO detection. By vary- 
ing t and r, a 2D-ES provides the frequency-controlled information 
of the first and fourth pulses. Hence, it is possible to make the first 
and fourth pulses sufficiently broadband that their frequency compo- 
nents at the ag and Pg transitions are of similar magnitude. Then, 
the sixteen ID experiments can be replaced with four 2D-ES where 
the second and third pulses are frequency controlled. A caveat in this 
identification is the assumption that the optical coherences evolve in 
a form like Eq.[TT] without errors of coherence transfers (SI-VII,IX). 

(c) Extension to overlapping pulses. The discussion above 
assumed negligible pulse overlaps. Remarkably, Eqs. ( 114b . l ll5b and 
( 116b still hold in general for any T,t > and T > 3a, with the 
exception that CJ^ in Eq.[T4]must be replaced by 



^ + ^'■^12^ + 2 



17] 



to account for the fact that the third pulse must act in the sample 
before the LO can detect the polarization (SI-IV-B). In the case of 
well-separated pulses, t ^ a, Eq.ll7lreduces to C^gCf,^. 

The measurements of the real and imaginary part of the PE sig- 
nal in the T,t = limit are recognized with the names of transient 
dichroism (TD) and transient birefringence (TB), respectively 1241 . 
and are very interesting for QPT. For resonant TD/TB (ug + tJ4 = 
ujrs + tt'sg), Eq. ( 117b reduces to ^CJ^^C^^, which shows that the 
LO monitors only half of the original polarization since it interferes 
with the polarization as it is generated. Consider such a resonant 
TD/TB experiment where, even though the pulses can achieve fre- 
quency selectivity, they are short in the sense that a j^, where A 
is a characteristic reorganization energy scale of the bath. In this sit- 
uation, the bath state will not evolve during the action of the first 
two pulses, allowing unambiguous preparation of initial excitonic 
states tensored with the same initial equilibrium bath configuration 
(SI- VIII), yielding a consistent QPT. Also, for T,t — 0, the free evo- 
lution of the optical coherences does not contribute to the signal, and 
the short timescale a does not allow for errors of population or coher- 
ence transfers to occur in the preparation or detection stages. Hence, 
a highlight of the TD/TB signal is that it is determined exclusively 
by the dynamics of the single-exciton manifold. Scenarios where the 
TD/TB configuration is preferred compared to the PE are excitonic 
systems coupled to highly non-Markovian baths 1251 1261 l27l 1281 . 

(d) Numerical stability of QPT. An investigation of the stabil- 
ity properties of the matrices associated with the reconstruction of 
x{T) from the sixteen enumerated experiments shows that our pro- 
tocol is very robust upon the variation of the structural parameters of 
the system, namely, the ratio between the two dipole norms ds /cIa, 
the angle between the site dipoles (/>, and the mixing angle 6. Gen- 
eral exceptions occur at the vicinity of 6 = 0, ^ where the coupling 
J vanishes, as well as for = and ds/dA — 1, that is, the 

homodimer case (SI-XI-A and B). 

Numerical example To test the extraction of x from experimental 
spectra, we consider a dimer with Hamiltonian parameters that are on 
the order of previously reported experiments consisting of light har- 
vesting systems |29 18] (cja = 12881 cm~\ ojb = 12719 cm"\ 
J = 120 cm~^, yielding 9 — 0.49). We assume a toolbox of two 
carrier frequencies uj+ — 13480 cm^^ and uj- — 12130 cm^^. 



-'By comparing Eqs. |7] and [14] we notice that botti lSpE]%i!s2-'<'3,'^i ^'^•'^'^f ^"'^ 
^ei'^e^^e^ ^""6 related to iPpE via a real proportionality constant. Alttiougti 

we stiali denote P^'^'^^,^^^ (t, T, t) loosely as a polarization, when referring to its real 
and imaginary parts, we must remember that they are proportional to the real and imaginary 
parts of the signal [SpE]^\]e2'e^,e'^ ('^^ ^rid to the imaginary and real parts of the 
actual polarization P pe^ respectively. 
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respectively, so that cji £ {ci;+,lj_} for all i, and the width of the 
pulses to be FW HM = 28.3 fs in intensity, which corresponds to 
(7 = 40 fs in amplitude. The parameters satisfy the MDC condition 
with C'/C" = 20. The pulses are long enough to guarantee the se- 
lectivity of the produced exciton, but short enough to allow for the 
evolution of the bath induced excitonic dynamics to be monitored. 
We choose ds/dA = 2 and (j) — 0.3. We present simulations on 
the QPT for this system, where each chromophore is linearly cou- 
pled to an independent Markovian bath of harmonic oscillators. The 
dissipative effects are modeled through a secular Redfield model at 
temperature T = 273 K (SI-X). 

Since we are working in the MDC regime, the signals in Eq. [14] 
are simply proportional to (-Pe'^e'atea.ei (t, T, t))i3o. Fig. 3 plots 
the sixteen real and imaginary parts of the {P5!z'z',z(0, T, 0))iso val- 
ues, which can be regarded as signals from the TD/TB setting, or as 
PE signals with the coherence and echo time propagators factored 
out. They have been calculated via an isotropic average of Eqs. il5l 
and J16b (SI-VI,XI-B). In our simulations, we consider an inhomo- 
geneously broadened ensemble of 10,000 dimers with diagonal dis- 
order. The site energies are drawn from Gaussian distributions cen- 
tered about (jJA and ujb, respectively, both with standard deviation of 
o'inh ~ 40 cm^^. For every waiting time T, the signal is calculated 
with this fixed ensemble. After a normalization step, the signals are 
of 0(1) or smaller. Additional noise simulating experimental errors 
due to laser fluctuations is included. This consists of independent re- 
alizations at every waiting time T of Gaussian noise on the measured 
signals with zero mean and aiaser = 0.05 standard deviation. 

Fig. 3 plots the ideal and inhomogeneously broadened, noisy sig- 
nals as continuous and discrete points, respectively. The ideal signals 
are calculated from a single dimer with no disorder and without laser 
fluctuations. All plots start at T = 3o- = 51 fs, since for earlier 
times, the initial states are not yet effectively prepared. Errors in 
experimental signals {Pz]i'Jz'"z{0, T, Q))iso translate into errors of re- 
constructed xiT). An estimate of the amplification of relative errors 
is set by the condition numbers of the matrices to be inverted, which 
is lower than k = 2.9 for our set of parameters ds/dA, 0, (j> (SI-XI- 
B). Reconstruction of x must respect the known symmetries, Eqs. 
I2I4I Eqs. l2l3l are built into the corresponding matrix equations, but 
Eq. must be included by using a semidefinite programming rou- 
tine. The latter is implemented using the open-source package CVX 
1301 , and the result is Fig. 4, which shows the discrete points repre- 
senting the reconstructed elements of x{T) from noisy data on top 
of the ideal results plotted as continuous functions of T. The relative 
error of the inverted x{T) averages to 0.12. Notice that despite the 
significant inhomogeneous broadening and noise, there is remarkable 
agreement between the ideal and the reconstructed values. This find- 
ing is reminiscent of studies due to Humble and Cina | 31 1. 

Fig. 4 illustrates the final objective of a coupled dimer QPT, 
namely, the process matrix x{T). Each panel shows processes of 
the density matrix as a function of T, conditional on the initial state 
being eta (a), (b), or afi (c and d), with the ideal T-dependence. 
The detailed balance condition in the Redfield model implies that 

XPfiaa = Xaappe.'^"'^^''^—, whlch cau be seen in panels (a) and 
(b). Also, note that due to the secular approximation, coherence-to- 
population and the reverse processes are zero. Note that all the deco- 
herence processes in our model occur within a timescale of hundreds 
of femtoseconds, with the afi coherence evolving through about three 
periods before practically vanishing (c and d). 

Clearly, a more complex interplay of the excitonic system with 
the phonon bath is possible ||321, but this example illustrates the 
essence of the type of information which can be obtained through 
QPT. 

Conclusions 

In this article, we have introduced QPT as a powerful tool to system- 
atically characterize the dynamics of excitonic systems in condensed 
phases. We identified the coherence, waiting, and echo intervals of 

the PE experiment with the state preparation, free evolution, and de- 
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tection stages of a QPT. In order to achieve selective state preparation 
and detection, we suggested frequency control through pulses of two 
different colors, although scenarios with time delays and pulse po- 
larizations as control knobs were also discussed here and elsewhere 
1331 1201 . By choosing between these colors for each of the four 
pulses, sixteen experiments can be carried out, which yield all the 
elements of x{T) related to the single-exciton manifold. An analy- 
sis of the reconstruction of x{T) in the presence of inhomogeneous 
broadening and experimental noise was provided, and the simulation 
on a model system shows that QPT of an excitonic system in con- 
densed phase is a very plausible goal. 

Equipped with x{T), which completely characterizes the exci- 
tonic dynamics, a plethora of questions can be rigorously answered 
about it. Some examples are: Can the bath be described as Marko- 
vian? If so, does the secular approximation hold, or can a population 
spontaneously be transferred to a coherence? t34J If not, what is its 
degree of non-Markovianity? 1351 1361 Does a given master equa- 
tion accurately describe the dynamics of the system? What is the 
timescale of each decoherence process? Are the baths coupled to 
each chromophore correlated? f37"381 How much entanglement is 
induced in the system upon photoexcitation? 1 39 1 Once these ques- 
tions are answered, interesting questions of control [40 | and manipu- 
lation of excitons can be asked. 

In summary, a QIP approach to nonlinear spectroscopy via QPT 
offers novel insights on the ways to design experiments in order to 
extract information about the quantum state of the energy transfer 
system. We believe this work bridges a gap between theoretical and 
experimental studies on excitation energy transfer from the QIP and 
physical chemistry communities, respectively. 
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Figure 3. Polarization signals from sixteen two-color experiments for fixed 
T = t = 0. The legends pqrs correspond to the real parts of the isotropically 
averaged signals (f ?;I',z',z{0, T, O))iso. The panels are organized by OPT 
initial state: (a) for |q)H, (b)for |/3)(/3|, (c) for |a)(/3|, and (d) for |/3)(q|). The 
ideal signals are depicted as continuous functions whereas the simulations with 
inhomogeneous broadening and noise are represented as discrete points of the 
same color. 
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